Localization of shadow poles by complex scaling 
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Through numerical examples we show that the complex scaling method is suited to explore the 
pole structure in multichannel scattering problems. All poles lying on the multisheeted Riemann 
energy surface, including shadow poles, can be revealed and the Riemann sheets on which they reside 
can be identified. 
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Resonant states are solutions of the Schrodinger equation with outgoing asymptotic boundary condition. It was 
pointed out long ago that these solutions must belong to complex eigenenergies, and the scattering matrix has poles 
at these energies jj]. In the coordinate space resonant eigenfunctions show oscillatory behavior in the asymptotic 
region, with exponentially growing amplitude, ~ exp[i(n — ij)r] (k, 7 > 0), thus they are not elements of the L 2 
I/"") \ space. Complex scaling is a most powerful and easily applicable method to describe such states It has been 

successfully applied in atomic |||| and nuclear physics . 

In single-channel problems the working mechanism of the complex scaling method (CSM) is well understood and 
there is almost no obscure point. The situation is not so clear, however, in multichannel cases. The pole structure 
of a multichannel scattering matrix is much more complicated than that of a single-channel 5-function. In the case 
of Hcrmitcan potentials a pole that would appear in one of the N channels in a single-channel problem, gives rise 
ly-j • to 2 7V_1 poles on different Riemann energy sheets in the coupled TV-channel problem ||||. The easiest way to label 
a Riemann sheet is to give the signs of the imaginary parts of the channel wave numbers ki (i = 1, 2, . . . , N) in an 
.TV-term sign-string (sign(Im k±), sign(Im /«2),...,sign(Im fcjy)) fl(|- In the zero-coupling limit all 2 N ~ 1 poles are at the 
same energy position on different sheets, while, by varying the coupling strengths, the poles move, and a crossing of 
the real energy axis by one of them implies a crossing over to another Riemann sheet. It has been a long-standing 
belief that only those poles (named ordinary poles) can have appreciable effects on the physically observable quantities 
which are on the Riemann sheet adjacent to the physical one. Recently, the effect of other poles, the so-called shadow 
pol es, of the multichannel scattering matrix on some physical observables has attracted interest in atomic pjj j, particle 
[[121], and nuclear physics [lC|Jl^j . It turned out that in certain cases the shadow poles can cause strong effects. E.g., it 
is a shadow pole which causes the very large cross section of the famous d + t—*a. + n thermonuclear reaction_ 
The effect of a shadow pole on the scattering matrix depends crucially on which Riemann sheet it is situated 

Although most of the applications of the CSM are in multichannel problems, up till now all investigations have 
been concerned with ordinary poles. In this Brief Report we show that, using the CSM, one can search for poles on 
5^ , different Riemann sheets and can identify the poles by their sheets. 

In a one-channel case the essence of the CSM is as follows. Instead of the 

H\V) = (T + V)\V) = E\V) (1) 
eigenequation of the Hamiltonian H, we solve the eigenvalue problem of the transformed Hamiltonian Hg = 

u{0)Hu- 1 {ey. 

H e \*e) =E g \y g } (2) 

(the 9 subscript of \& means that the wave function implicitly depends on 9; Hamiltonians with different 9 results in 
different wave functions). U(9) is an unbounded similarity transformation [Q, which, in the coordinate space, acts 
on a function f(r) such that 

U(6)f(r) = e™' 2 f{re* e ). (3) 

(If 9 is real, U{9) means a rotation into the complex coordinate plane, if it is complex, it means a rotation and scaling.) 
The two problems are connected by the ABC theorem fsi. If V is a (dilation) analytic operator: 
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(i) The bound eigenstates of H are the eigenstates of Hg, regardless of the actual value of 9, within < 9 < tt/2; 

(ii) the continuous spectrum of H will be rotated by an angle 29; 

(iii) a complex generalized eigenvalue of Eq. (^), E IeB = e — i\T, e,T > (with the wave number k TCS = n — vy, 
k, 7 > 0), belongs to the proper spectrum of Hg provided 28 > | axgE ies \. 

Roughly speaking, the complex scaling transformation changes the asymptotic wave function from exp[i(n — ij)r] to 
exp[i(« — vy)r exp(i6)], which, in the case of 29 > | arg-E res | = 2| argfc res |, makes the diverging wave function localized. 
It is important to note that, if the sign of fc rcs is reversed, then the outgoing wave with 9 = is localized, and the 
complex scaling spoils the localization unless 29 < | arg E ICS | . In a single-channel problem with a Hermitean potential 
— fc res is on the physical sheet, where there are no resonance poles, but in a multichannel problem the imaginary parts 
of some channel wave numbers may be positive. 

In multichannel cases Eq. ([I]) becomes a matrix equation 

JV 

Y J H a p\^)=E a \^ a ), a =1,2,..., AT, (4) 
0=1 

where the Greek letters are the channel indices. The transformation operator of Eq. (|^) becomes 

U aP (6) = S aP U(9). (5) 

In the literature there are some hints on the strange behavior of the CSM in multichannel cases. E.g. in jl6j the 
authors found that varying the rotation angle 9, channel continua can absorb resonances that were revealed before, 
however, they did not explain this phenomenon. In Q it is stated that such a phenomenon becomes transparent if 
one studies the multichannel problem on the Riemann energy surface, but no attempt has been made to assign these 
poles to Riemann sheets. 

Here we study the working mechanism of the multichannel CSM in a simple model, which is easy to comprehend 
and control. Our model consists of a target with two internal states, whose thresholds are E\ and E2, respectively, 
and a projectile. The two target states are the two channels. We choose one-term separable potentials |l7]] for both 
the diagonal and interchannel interactions, 

V a0 = \<p%(b))\ a p(4(b)\, a,/3 = l,2, (6) 



where \ipo(b)) is the eigenfunction of the three dimensional harmonic oscillator with n = oscillator quantum, b is 



the size parameter and A Q /3 are the (real) potential strengths (A12 = A21). For the wave functions \&g) of (0), we use 



the following trial functions 



I^Hj^ktf®), 7 = 1,2 (7) 



i=0 



in a variational method for the expansion coefficients Ci (this is the well-known wave function expansion (WFE) 
method). With this ansatz, functions Y$>1) 6 L 2 are selected. The use of harmonic oscillator functions both in (||) 
and ([7]) makes it possible to calculate all necessary matrix elements analytically |l7j]. We choose size parameters in 
(0) different from that in (0) so as to make the trial function more flexible. If we set the strength of one of the 
diagonal interactions to be zero, the selected channel cannot accommodate a resonance, so that all poles we find in 
the coupled-channel problem must originate from the other channel, which implies that the sign belonging to the 
other channel must be negative. For the sake of simplicity, we take s-wave states throughout; although, the analytical 
expressions of JlTj can be used for I ^ as well. 

As a first example, we choose A22 = 0, An = 1.0, b = 0.6, b = 2.0, E\ = 0, and E2 = 2 (we use atomic 
units H = m = 1). The basis sizes (ni, 712) are chosen so as to reach stable convergence. In the uncoupled case 
(A12 = A21 = 0) only one resonance pole appears at E = 3.049 — i2.153 [Fig. |l|(a)]. Switching on the coupling 
(A12 = A21 = 1.0), we get different pole arrangements at different rotation angles [Figs. |l](b)-(d)]. We can see in all 
figures that the poles distort the continua, as if they attracted or repelled the continuum points. This phenomenon 
has surfaced several times earlier, e.g. Refs. [p|]I§|], but as far as the author knows, it is yet unexplained. Furthermore, 
we emphasize that in this work the central question is the working mechanism of the CSM in multichannel problems, 
which requires the use of very different rotation angles. Thus we do not perform an optimization in the 9 parameter 
(which could be done by choosing the stationary point of the ^-trajectory, see e.g. M). These figures show the 
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appearance and disappearance of poles, the same phenomenon as was mentioned above. The choice A22 = and the 
fact that in the uncoupled case there is only one resonance pole guarantees that in this problem there is a pole on 

the ( ) Riemann sheet and another, a shadow pole, on ( — h) ||. The pole at 4.742 — «1.810 is revealed only when 

both continua have swept over this point, which implies that the condition 29 > | axgE^ es \, 7 = 1,2 must be fulfilled 

in both channels (where E^ es and E^ es are the channel energies). Consequently, this pole is on the ( ) sheet. For 

the other pole, which is on the ( — h) sheet at 2.395 — HA67, the relations 29 > | &rgE} cs \ and 29 < \ &igE^ cs \ must 
hold (cf. the remark after the ABC theorem, above). This is in full agreement with what we can see in the figures. 

This example tells us that, if a pole is revealed, then it is on a Riemann sheet which is characterized by negative 
signs for all channels whose continua have swept over the pole, and positive signs for all others. A pole at E is an 
ordinary one (i.e. it is on a sheet adjacent to the physical sheet) if it has been swept over by the continua of all 
channels whose threshold energies are lower than He(E) and has not been swept over by any other ones. From this 
it follows that one can imagine situations where a shadow pole can be revealed only if the rotation angles in different 
channels are different. For instance, to reveale a pole above the first channel threshold on the (+— ) sheet, the rotation 
angle in the second channel must be greater than the one in the first channel. It is questionable whether the CSM can 
cope with such a constraint. The theory of the multichannel CSM always assumes that 9 is the same for all channels. 
Here, just to see what happens, we venture to choose two different 9. 

It seems natural to generalize the multichannel complex scaling transformation (|^) in the following way: 

U a0 =6 af3 U a (9 a ). (8) 

In the coordinate space the action of U a (9 a ) is 

U a (e a )f(r) = e 3ie «/ 2 f(re ie «). (9) 

This definition ensures that the {U a p\ operator matrix inherits all properties of U{9). Applying the transformation 
(0) to the Hamiltonian of Eq. (Q) , we arrive at 

N 

Y J U a {e a )H aP Up 1 {9p)\¥^) = 

13=1 

El >e0 \^J, a = 1,2,..., N. (10) 

As an illustrative example, we write down the function U a (9 a )H a pUp 1 {9p)\^g ii ) and its overlap with the function 

(<Fg I in the coordinate space. If the operator H a p connects channels which have the same dynamical coordinate r, 
then 

(r\U a {9 a )H aP Up 1 {9p)\^ 9f ) = 

e 3i(9--flfl)/2 1 r ( ^( re tf«.)^( re i(» B -^)) j (U) 

and 

x / njre-^H^r)^ (r e -^)r 2 dr. (12) 



If the H a fj operator connects rearrangement channels with the dynamical coordinates r a and rp [i.e. in the coordinate 
space H a pf(rp) = J drpr 2 H a fi(r a ,rp)f(r^)], then 

{r a \U a {0 a )H aP U^{9 p )\^) = e ^-M/2 

x J H a0 (r a e ie %r^ gfi {r e- l ^)rldrp, (13) 

and 

{^U a (9 a )H aP Up\9p)M) = 



xV? (rpe-^yldr^ldrp. (14) 
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Deriving ( |l2| ) and (|lj) the Cauchy theorem was used, assuming that all potential operators are analytic [like ours (0)]. 
We can see that the matrix elements of the transformed operator U a {9 a )H a ijU^ 1 {9p) between the original channel 
states |*£ o ) and ) can be expressed as the matrix elements of the original operator H a p between the so-called 
back-rotated channel states. This is a well-known feature of the usual complex scaling method, too. 

We tested this generalized CSM with the above two-channel problem setting 9 2 = 0. The result is in Fig. ^. The 
position of the ( — h) shadow pole remains the same as was in Figs, [j] within 7 decimal digits, which is a remarkable 
stability regarding that no optimization was made in b and 9±. 

The really relevant test is, however, an example where there is a shadow pole on (-1 — ). To achieve this, we set 
A22 = 0.2, which, in a one-channel problem, results in a pole at 2.430 — i3.704 [Fig. ||(a)]. Switching on the coupling 
(A12 = A 2 i = 0.42), figures similar to Fig. 1 can be generated [Figs. |](b)-(d)]. Now the pole at 3.381 — z3.228 is 

revealed when both channel continua have swept over this point, so that this pole is on the ( ) sheet. The other 

pole at 2.373 — 13.357 is revealed when the continuum of the second channel has swept over it and that of the first one 
has not, which shows that this pole is a shadow pole on the (H — ) sheet, in agreement with the fact that these poles 
originate from the second channel H . In this example the variation of the rotation angles slightly removes the poles 
from their original positions. This is, however, certainly caused by the fact that, because of the unlucky location of 
the poles, we have to choose rotation angles that are far from optimum. If we do not want to reveale the two poles 
at the same time, we can optimize the 9 angles, which results in stable pole positions. 

Finally, we mention an interesting feature of the present method. Let us suppose that there is a multichannel problem 
where there are degenerate thresholds. Then some of the Riemann sheets cannot be reached from the physical sheet 
by following analytical continuation paths because one cannot pass between two thresholds that coincide. Using the 
above multichannel CSM with different rotation angles in these channels, we can reach such Riemann sheets. 

In summary, we have investigated the applicability of the multichannel complex scaling method to explore the 
pole structure in multichannel scattering problems. We have used a natural extension of the single-channel complex 
scaling transformation to multichannel cases, which allows us to find the poles of the scattering matrix on all Riemann 
sheets. We have found that this extension works as expected and is able to find all conventional poles and shadow 
poles reliably. 
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FIG. 1. Energy eigenvalues of (a) a one-channel problem with a one-term separable potential (An = 1.0); (b)-(d) a 
two-channel problem with one-term separable potentials (An = 1.0, A22 = 0.0, A12 = A21 = 1.0, and E2—2). The dots 
are the points of the rotated discretized continua, while the circles are the poles of the S-matrix on different Riemann sheets. 
The rotation angles (in radians) are: (a) 0.4, (b) 0.2, (c) 0.4, and (d) 0.7. 



FIG. 2. Energy eigenvalues of the two-channel problem of Fig. hj. The rotation angles are: 9\ = 0.45 and 62 = 0. 



FIG. 3. The same as Fig. with the potential strengths: (a) A22 = 0.2; (b)-(d) An = 0.0, A 22 = 0.2, and A12 = A21 = 0.42. 
The rotation angles are: (a) 0.78, (b) 0.68, (c) 0.78, and (d) X = 0.42 , 9 2 = 0.78. 
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